Beyond-mean-field study of the possible "bubble" structure of '^'^Si 
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Recent self-consistent mean-field calculations predict a substantial depletion of the proton density 
in the interior of ^*Si. In the present study, we investigate how correlations beyond the mean 
field modify this finding. The framework of the calculation is a particle-number and angular- 
momentum projected Generator Coordinate Method based on Hartree-Fock-Bogoliubov-|-Lipkin- 
Nogami states with axial quadrupole deformation. The parametrization SLy4 of the Skyrme energy 
density functional is used together with a density-dependent pairing energy functional. For the first 
time, the generator coordinate method is applied to the calculation of charge and transition densities. 
The impact of pairing correlations, symmetry restorations and shape mixing on the density profile 
is analyzed step by step. All these effects significantly alter the radial density profile, and tend to 
bring it closer to a Fermi-type density distribution. 
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I. INTRODUCTION 



Charge distributions in atomic nuclei [l|-|7[ provide 
very detailed information about nuclear structure. They 
are obtained through the analysis of elastic electron- 
nucleus scattering data. Because of the absence of suit- 
able targets, data on unstable nuclei are available, to the 
best of our knowledge^^only for [| and [lll^. The 
SCRIT project |ll| - [l3| of constructing a high-resolution 
electro n sp ectrometer that is underway in Japan and 
ELISe planned to be constructed at FAIR, are ex- 
pected to provide data about the charge distributions 
and transition form factors for many exotic nuclei in the 
future. 

Because of the saturation properties of the nuclear 
medium, the radial dependence of the nuclear density 
takes, at the lowest order, the form of a Fermi distri- 
bution. However, the density often deviates from this 
simple behavior because of quantal effects related to the 
filling of single-particle states with wave functions that 
have specific spatial behavior. In this context, si/2 or- 
bits in spherical nuclei have a very peculiar signature, as 
they are the only ones that contribute to the density at 
the nuclear center. Depending on whether they are filled 
or empty, Si/2 orbits can generate a central bump in the 
density as it has been observed for ^°Ca [l5|, or a central 
depression. 

Mean-field-based methods [TB] are the tools of choice 
when modeling the nuclear density distribution. Indeed, 
they include the ingredients required for this task: the 
full model space of occupied single-particle states as de- 
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grees of freedom together with an effective interaction 
that reproduces the empirical saturation properties of 
nuclear matter. 

The density profile and the spatial dependence of the 
single-particle potentials are closely related and self- 
consistently linked to each other. A central depression 
in the density might be accompanied by two very spe- 
cific properties of the mean-field potential, one related to 
its central part and a second one related to the spin-orbit 
potential. 

A central depression of the density is reflected in the 
central potential by a maximum at the origin and a min- 
imum for some finite distance r. This is often called a 
"wine-bottle" shaped central potential, referring to the 
shape of the bottom of a bottle of wine. Levels with low 
orbital angular momentum ^ are then pushed up rela- 
tively to those with large i that are pulled down. For suf- 
ficiently large rearrangement, the order of single-particle 
levels can even change, lowering the central density even 
more and leading to the so-called "bubble nuclei". For 
specific "bubble magic numbers" 18, 34, 50, 58, 80, 120, 
. . . (itI . ITsI I large shell effects might compensate for the 
loss in binding energy due to the reduced central density 
well below the nuclear matter saturation value. There 
was some speculation in the 1970s whether such struc- 
ture could exist in nuclei that were about to become ac- 
cessible for detailed studies, in particular "^^Ar and some 
Hg isotopes [l7l - [23 |. However, the possibility of a bub- 
ble structure in these nuclei has been ruled out by ex- 
periment. By contrast, predictions that superheavy and 
hyperheavy nuclei beyond the currently known region of 
the mass tabic might take the form of bubbles [25l - [29| 
arc still standing. In fact, for a large charge number Z , a 
hollow density distribution is energetically favored over 
a regular one as it lowers the Coulomb repulsion. In this 
context, one often distinguishes between "true bubbles", 
which have vanishing density in their center, and "semi- 
bubbles" , which have a central density significantly lower 
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than saturation density, but with a non-zero value. 

The second effect of a central depression concerns the 
spin-orbit potential. In self-consistent mean-field mod- 
els, this potential is proportional to the gradient of a 
combination of proton and neutron densities, whose rel- 
ative weights depend on the model and parametrization 
[27} . For nuclei with a regular density profile, it is 
peaked at the nuclear surface. For nuclei with a central 
depletion of the density, the spin-orbit potential has a 
second peak of opposite sign in the nuclear interior. This 
usually reduces the spin-orbit splitting of orbits located 
mainly at the nucleus center, whereas that of orbits sit- 
uated at the nuclear surface is not affected. 

Recently, there has been a renewal of interest in nu- 
clei presenting a hollow in their density distributions. 
Some modern parametrizations of the relativistic mean 
field l30l. |3 l| an d of the Skyrme energy density functional 
(EDF) [32, 133 predict a hollow proton density for '^'^Si 
and some neutron-rich Ar isotopes. At the time being, 
•^^Si stands out as the only candidate on which many dif- 
ferent effective interactions agree. The possible proton 
bubble structure of this nucleus has also been suggested 
as an explanation of the results on the transfer reactions 
36S(d,p)37S and ^^Si{d,pf^Si. Indeed, the splitting be- 
tween the observed neutron 3/2^ and 1/2^ levels that 
have the largest spectroscopic factors in the 2p shell is 
decreased from 1.7 MeV) to ^^Si (w 1.1 MeV) 

mil. 

Besides the debatable interaction dependence of the 
density distributions, one may wonder whether bubble- 
type structures are stable against correlation effects. In- 
deed, any correlation will inevitably populate empty lev- 
els and in particular the 2si/2, even in models like the 
one that we use here, where its population cannot be 
easily singled out. In fact, it is known for a long time 
that, already for nuclei with a more regular density dis- 
tribution, mean-field calculations tend to overestimate 
the spatial fluctuations of the density when compared to 
data. Correlations usually tend to flatten out the den- 
sity distribution, and often bring it closer to data. The 
effect of pairing has been studied in Ref. [sB], and the 
impact of fiuctuations in shape degrees of freedom has 
been studied within the random phase approximation 
(RPA) for many spherical nuclei p37H44i a nd using a one- 
dimensional [il] or five-dimensional j39l. |40| . l46j micro- 
scopic collective Hamiltonian for some transitional ones. 

The most obvious correlations that could reduce the 
central depression of the density are due to pairing [2l|. 
However, many calculations made for bubble nuclei do 
not include them pj,, 20,, 22, ill] . This can be justified for 
'^'^Si, where the large gap between the proton 1^5/2 and 
2si/2 levels suppresses pairing, resulting in its unphysical 
collapse when treated with the Hartrce-Fock-Bogoliubov 
(HFB) method. Pairing correlations have then to be 
treated beyond the mean field approximation. Another 
kind of correlations that might affect the density profile is 
related to the spreading of the ground-state wave func- 
tion around the mean-field configuration. It has been 



pointed out in Ref. that the ground states of most 
light nuclei may show strong shape fluctuations that in 
general lead to a substantial increase of their charge radii 
when treated in a beyond-mean-field framework. The 
same effect may also strongly influence the density pro- 
file. 

In the following, we will compare results obtained 
from calculations that successively add correlations to 
the ground-state wave function: 

(i) spherical mean-field calculations without taking 
pairing correlations into account (HF) 

(ii) spherical mean-field calculations including pair- 
ing correlations within the HFB-|-Lipkin-Nogami 
(HFB-I-LN) scheme, which constitutes an approxi- 
mate variation after projection on particle number 

(iii) particle-number projection after variation of the 
spherical mean-field state obtained in HFB+LN 

(iv) configuration mixing of angular-momentum J = 
and particle-number projected mean-field states 
with different intrinsic axial quadrupolc moment. 
We will refer to these wave functions in the fol- 
lowing as symmetry-restored generator coordinate 
method (GCM). 

In addition, we will study how much the observable 
charge density, obtained through convolution of the pro- 
ton density with the proton's internal charge distribution, 
differs from the point proton density used to calculate the 
energy. 

Symmetry-restored GCM has been used to describe 
ground-state correlations and collective excitation spec- 
tra of a large range of nuclei with reasonable success 
[48l - l53| . Its actual implementations are not yet flexible 
enough to reproduce correctly all details of excitation 
spectra, mainly because of the usually too low moment 
of inertia. However, this method describes rather well 
properties related to the nuclear shape, such as transi- 
tion probabilities. It is therefore important to enlarge 
its range of applications by the calculation of charge and 
transition densities in the laboratory frame. In this pa- 
per, we report on a first application that enables us to 
illustrate the effect of various kinds of correlations on the 
density distribution of ^*Si. 

We will first give a brief outline of the model. Re- 
sults for low-lying collective states in "^^Si are presented 
in Sect. IIII[ whereas the modification of the ground state 
density distribution brought by correlations is discussed 
in Sect. ITVl Section |V] will summarize our findings. 



II. CALCULATIONAL DETAILS 

The self-consistent HFB equations arc solved on a cu- 
bic three-dimensional coordinate-space mesh extending 
from —11.2 fm to 11.2 fm in each direction with a step 
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size of 0.8 fm. Thanks to the reflection symmetry with re- 
spect to the a; = 0, y = and z = planes imposed on the 
single-particle wave functions in our code [5j| , it is suffi- 
cient to solve the HFB equations in 1/8 of the box. The 
HFB equations are complemented by the Lipkin-Nogami 
prescription to avoid the unphysical breakdown of pairing 
correlations at low level density. A constraint on the ax- 
ial mass quadrupole moment q = {Q2) = (Zz^ — — y^) 
is used to construct mean- field states \q) with different 
intrinsic deformation. 

Eigenstates of the particle-number operators N and Z 
are obtained by applying the projection operator 
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for neutrons and protons. Eigenstates of the total an- 
gular momentum in the laboratory frame and its z 
component with eigenvalues h'^J{J + 1) and HM, re- 
spectively, are obtained by applying the operator 



^MK ~ 



2J + 1 

87r2 



da 



dp sin(^) 
^0 



iiK R (2) 



that contains the rotation operator R 



-iajj: 



V g-n-h Wigncr rotation ma- 

trix 2?j(^^(a, /3, 7) on the nucleus' wave function. Both 
depend on the Euler angles a, /3 and 7. The operator 
^iiK picks the component with projection K along the 
intrinsic z axis from the mean- field state. Throughout 
this study, we will restrict ourselves to axial states with 
K = 0. As angular-momentum projected states will be 
always projected also on particle-number, we drop the 
indices iVo and Zofor the sake of notation: 



\JMq) = 



[q\pJoPNoPzo\q) 



(3) 



GCM 155| is a very flexible tool that in particular allows 
us to study the spreading of the mean-field ground-state 
wave function in collective degrees of freedom. It will 
be used here to study the fluctuations of the spherical 
ground state of '^'*Si with respect to the axial quadrupole 
moment assuming a superposition of projected HFB+LN 
states of different deformation \q): 



\jMp) = Y,r^M)\JMq) 



(4) 



The weight factors /^f (g) and the energies E'^ of the 
states I JMji) are the solutions of the Hill- Wheeler- Griffin 
equations j55| 



energy density functional that depends on the mixed den- 
sity matrix between the two projected states | JMq) and 
\JMq') H^. 

Throughout this article, we will use the parametriza- 
tion SLy4 |57| of the Skyrme energy density functional 
together with a local pairing energy functional of sur- 
face type with parameters po = 0.16 fm"'^ for the 
switching density and Vq = —1000.0 MeV fm'^ for the 
pairing strength unless noted otherwise. A soft cutoff at 
±5 MeV around the Fermi energy is used when solving 
the HFB equations as described in Ref. [5^ . More details 
about the calculations of the GCM kernels can be found 
in Ref. [sij and references given therein. 

The weight functions /^(f?) in Eq. ^ are not orthogo- 
nal. A set of ortho normal collective wave functions g^{q) 
can be constructed as [551 



Y.{Af-'f\q,q')f^{q'). 



(6) 



It has to be stressed, however, that does not 

represent the probability to find the deformation q in 
the GCM state \JMiJ.). In addition, in the absence of 
a metric in the definition of the correlated state | JM/i) , 
Eq. (HI, the values of g^(q) for a converged GCM solution 
still depend on the discretization chosen for the collective 
variable q, which is not the case for observables like the 
energy or transition probabilities. 

The spatial density distribution of the projected GCM 
states is constructed as the expectation value of the op- 
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87r2 
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X {q\R[P) E n) PnoPzo W) , 



(7) 



where we use the shorthand = (a, /?, 7) for the Euler 
angles. Note that the calculation of the density in the lab- 
oratory frame requires projectors on the left and on the 
right. More details on the calculation of the correlated 
density will be given in a forthcoming publication [59^ . 



J2[n'iq,q')-Elf^-'iq,q')] //(g') = , (5) 



III. SPECTROSCOPY OF LOW-LYING STATES 



where the norm kernel reads N''{q,q') = {JMq\JMq') 
and where the energy kernel is given by a multi-reference 



The energy curves obtained after projection on par- 
ticle number and on angular momentum are shown in 
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FIG. 1: (Color online) Energy curves for the particle- 
number-projected HFB states (N&Z) and particle-number 
and angular-momentum projected states (J = 0, 2, 4, 6 
curves) for "^"'Si as a function of the intrinsic quadrupole defor- 
mation of the mean-field states they are projected from. The 
solid square dots correspond to the lowest GCM solutions, 
which are plotted at their average deformation ^ q \g'^{q)\'^ 
(see text). 
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FIG. 2: (Color online) Nilsson diagram of the eigenvalues of 
the single-particle Hamiltonian for neutrons (left panel) and 
protons (right panel) as obtained with the Skyrme interaction 
SLy4 for ^'*Si as a function of the quadrupole deformation. 
Solid (dotted) lines represent levels of positive (negative) par- 
ity, and black, red, green and blue color represents levels with 
expectation values of (jz) = 1/2, 3/2, 5/2 and 7/2. The thick 
long-dashed line represents the Fermi energy. Single-particle 
levels are labeled for the spherical configuration only. 



Fig. [TJ The abscissa corresponds to the mass quadrupole 
moment q of the intrinsic states that is projected (upper 
scale) and, equivalently, to the dimensionless quadrupole 
deformation 

where R = 1.2 A^/'^ fm. 

The particle-number-projected energy curve presents 
a spherical minimum with a steep rise with deformation 
(dotted line), as expected for a nucleus with large neu- 
tron and proton shell gaps, cf. Fig. [5] The projection on 
total angular momentum J = leads to energy curves 
with prolate and oblate minima at about the same defor- 
mation 1/32 1 ~ 0.26. The presence of these two minima 
is a usual result of angular-momentum projection when 
the no n-proje cted energy curve presents a spherical min- 
imum j47H49| . In fact, at small deformation, the states 
of a given J projected out from prolate and oblate mean- 
field states with the same |/32| value are almost equiva- 
lent. Oblate and prolate minima are also obtained for 
higher J-values. Our results are very similar to those of 
a similar calculation using the Gogny interaction [60| . 

The energies E'^^ of the lowest GCM states are also 
shown in Fig. [T] (solid square with a label J^) at the 

mean deformation \g'l{l)\^ of the mean-field states 
on which they are build. This mean deformation is 
not an observable; still, it often provides a good indica- 
tion about the dominating mean-field configurations in 
a GCM state. The mean deformations and the B{E2) 



transition strengths suggest to organize the correlated 
states into the two structures displayed in Fig. [Sj where 
the computed transition probabilities and the energy of 
the levels are also compared with the available experi- 
mental values [6l[. Our result can be interpreted as re- 
sulting from the coexistence of an anharmonic spherical 
vibrator and a prolate deformed band at low excitation 
energy. Both structures are not pure and distorted by 
their strong mixing. 

The energy of the recently observed low-energy Oj' 




Exp. SLy4 



FIG. 3: (Color online) Comparison between the experimental 
(left) and calculated (right) excitation energies Ex and B{E2) 
values (in units of fm^) for the low-lying states of ^*Si. 
Experimental data are taken from Ref. [gH . 
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FIG. 4: (Color online) Collective wave functions (cf. 
Eq. de))) of the two lowest GCM solutions for J = 0. 



state pT| and the out-of-band B{E2 : 2j'" — > Oj) value are 
reproduced rather well. However, the electric monopole 
p^iEO]0+ 0+) and the in-band B{E2;2+ 0+) are 
overestimated by our model: 58.1 x 10~^ compared to the 
experimental value of 13.0 (0.9) x 10~^ [6l| for the former 
and 60.5 fm^ compared to 16.6 fm'' for the latter. 
This discrepancy might indicate [63 that the two lowest 
0+ GCM states are too strongly mixed in our calcula- 
tion. The corresponding collective wave functions g^{q) 
are displayed in Fig. |4l Both are indeed spread over a 
very wide range of deformations, with similar contribu- 
tions at small deformation |/32| ~ 0. The ground state is 
peaked around the deformations of the two minima in the 
J = projected energy curve, cf. Fig.[T] By contrast, the 
wave function of the Oj state is peaked at large prolate 
and oblate deformations where at least one downsloping 
level from the neutron shell becomes intruder by 
crossing the upsloping levels from the sd shell, cf. Fig. [2] 
This is consistent with the interpretation of the Oj state 
in '^^Si as a counterpart of the deformed ground state of 
the slightly lighter nuclei located in the so-called "island 
of inversion" [61| . 



IV. DENSITY DISTRIBUTION 

To quantify the depletion of the proton density distri- 
bution, we will use a depletion factor 



Pmax.p Pccnt.p 
PmsiX^p 



(9) 



which has been used in Ref. [3l|, |33| and that measures 
the reduction of the density at the nucleus center rela- 
tively to its maximum value. 

The effect of pairing correlations, projection on good 
quantum numbers and configuration mixing on the ra- 
dial profiles of the proton, neutron and total densities is 
displayed in Fig. [S] The densities of the HF, HFB-hLN 
and particle-number projected HFB-I-LN states are com- 
pared to those of the GCM O"*" ground state. To facil- 
itate the comparison, the proton and neutron densities 
are rescalcd hy A/Z and A/N factors, respectively. 



A large depiction at r = and a bulge at r 1.8 fm 
are obtained for the proton density when the HF method 
is used (top left panel of Fig.[5|). The HF neutron density, 
however, has an opposite behavior, with a flat shoulder at 
intermediate r values and a bump at the nucleus center. 
This bump is similar to the one found experimentally for 
the charge density in ^°Ca [l5|. Altogether, the total 
density has an almost flat, even slowly rising, profile in 
the interior of the nucleus. The same compensation of 
neutron and proton densities in the system's interior is 
also found at all other stages of the calculation. 

Unconstrained HFB calculations for ■^'*Si give the same 
result as the HF approximation. This is because the large 
Z = 14 gap of about 4.5 MeV between the proton 1^5/2 
and 2si/2 levels in the single-particle spectrum prevents 
the protons from becoming superfiuid at the HFB ap- 
proximation. The even larger A'^ = 20 gap in the single- 
particle spectrum has the same effect for neutrons. The 
situation is different for nuclei such as ^^O and ^^Ar, 
where pairing correlations are already active at the HFB 
level and wash out the bubble structure predicted by HF 
calculations [33| . 

The collapse of pairing correlations when the density 
of single-particle levels falls below a critical value is a de- 
ficiency of the HFB method [H, [11] . It can be partially 
corrected by using the LN procedure (top right panel 
of Fig. O. The level occupation is then smeared over 
the Fermi energy and the proton 2si/2 orbital becomes 
partially occupied. As a consequence, the central pro- 
ton density rises considerably from F^ax = 0.41 (HF) 
to Fmax = 0.32 (HFB-I-LN). The HFB+LN density pre- 
sented in Fig. O is calculated using occupation numbers 
corrected for particle-number projection by the approx- 
imation described for example in Ref. |64l |. Using the 
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FIG. 5: (Color online) Neutron, proton and total radial den- 
sities at X — y — for ''''Si for the spherical HF state (top left 
panel), the spherical HFB+LN state (top right) and its pro- 
jection on good particle numbers {N&[Z, bottom left), as well 
as for the GCM 0^ ground state {N&cZ, J = 0, GCM, bottom 
right). Neutron and proton densities have been rescaled with 
the factors A/N and A/Z, respectively. 
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FIG. 6: Contour plots of the neutron (upper panels) and pro- 
ton (lower panels) densities in the y = plane for the par- 
ticle number projected HFB-I-LN state with /32 = 0.26 (left 
column), its projection on both particle numbers and total 
angular momentum J — (middle column) and for the 0"'' 
GCM ground state (right column). 



non-corrected BCS occupation numbers instead would 
overestimate the effect of pairing and give a mucli larger 
reduction of the depletion factor. 

Projection of the HFB-I-LN state on good particle num- 
bers (bottom left panel of Fig. [5|) substantially reduces 
the pairing correlations, and the density profiles almost 
go back to the HF ones with i^max = 0.36. This reflects 
the well-known fact that the LN approximation overesti- 
mates the correlations in the weak pairing limit (whereas 
HFB underestimates them) , cf . for example Ref . [6a] , and 
indicates that in this case a correct treatment of pairing 
requires to go beyond the mean field. 

The behavior of the density close to the origin is usu- 
ally discussed in terms of the occupation of single-particle 
states in the spherical HF basis. This is not obvious in a 
method like the one that we use where the mean-field ba- 
sis is different for each deformation. Deformation mixes 
single-particle states with different orbital angular mo- 
mentum. In particular, when one expands a deformed 
basis in terms of the spherical one, the proton 2si/2 level 
gets partially filled. The situation is even more compli- 
cated after projection and configuration mixing, cf. the 
bottom right panel of Fig. [5] We have seen in Fig. 0] that 
the collective wave function of the 0"*" GCM ground state 
is spread over a wide range of intrinsic deformations. 

Figure[6]illustrates how the density distribution of neu- 
trons (upper panels) and protons (lower panels) is modi- 
fied at different levels of our calculation. The left column 
shows contour plots of both densities for the particle- 
number projected HFB+LN state with f32 — 0.26 that af- 
ter angular-momentum projection gives the prolate mini- 
mum of the J = curve in Fig.[TJ The proton density still 
exhibits a central depletion, but less pronounced than it 
is for the spherical HFB-I-LN state, reflecting the partial 
fllling of the 2si/2 level by deformation. After projection 
on total angular momentum J = (middle column), the 



density is obtained in the laboratory frame and is spheri- 
cal. However, the central depression of the density is sim- 
ilar to the one found when projecting on particle numbers 
only. The configuration mixing leading to the GCM 0+ 
ground state (right column) increases the central proton 
density again and simultaneously reduces the value at 
the bulge, see the bottom right panel of Fig. \5\ which in 
combination reduces the depiction factor to -Fmax = 0-21. 
The values of central and maximum densities and of the 
depletion factor for the states discussed above are sum- 
marized in Table m 

Up to now we have discussed the density of point pro- 
tons and neutrons, thereby neglecting that protons and 
neutrons are composite particles of extended size. When 
searching for experimental signature of a depleted cen- 
tral density in '^^Si by elastic electron scattering, how- 
ever, this has to be taken into account. The observable 
charge density is calculated by convolution of the proton 
density with a Gaussian form factor ,66] with a proton 

size a = \/2/3 {r'^)']/^ = 0-65 fm, which for spherically 
symmetric density distributions leads to 



Pchir) = 



awn 



dr'r' Pp{r') 



-(r-r')Vi^ 



(10)' 

The charge density (right panel) is compared in Figure [7] 
to the point proton density (left panel) for the same four 
cases discussed in Fig. [5l Like correlations, the convolu- 
tion ([To)) tends to even out the variations of the density 
profile: the central density rises and the maximum den- 
sity of the outer bulge becomes smaller. In combination, 
both leads to a substantial reduction of the depiction fac- 
tor from i^inax = 0.41 for the point proton density in a 
spherical HF calculation to i^max = 0.09 for the charge 
density of the 0"'" GCM ground state. 

Adding correlations, the root-mcan-square (rms) ra- 
dius of the point proton density increases from 3.127 fm 




3 4 
r(fm) 



FIG. 7: (Color online) Comparison of point-proton densities 
(left) with the folded charge densities (right) for ^'*Si for the 
same states as in Fig. [S) 
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TABLE I: Central and maximum proton density of Si 
and the depletion factors -Fmax.p [cf. Eq. Q] and Fsat.r [cf. 
Eq. For the latter, values for proton, neutron and the 

total densities are given. All densities are in fm"^. The three 
first lines correspond to a spherical state. The values labeled 
with N&iZ, J — correspond to the prolate minimum of the 
N&iZ, J — projected energy curve of Fig. [T] 







pmax,p 


max,p 


-^sat , p 




-^sat ,£ 


HF 


0.044 


0.074 


0.41 


0.34 


-0.37 


-0.08 


HFB+LN 


0.050 


0.074 


0.32 


0.24 


-0.31 


-0.08 


NSzZ 


0.047 


0.074 


0.36 


0.28 


-0.30 


-0.06 


N&iZ, J = 


0.051 


0.073 


0.30 


0.22 


-0.27 


-0.07 


GCM(g.s.) 


0.057 


0.073 


0.21 


0.13 


-0.22 


-0.07 



for the spherical HF state to 3.133 fm for the particle- 
number projected spherical HFB+LN state and to 
3.180 fm for the GCM 0+ ground state. Looking at the 
density profiles in Fig. [7] the larger radius of the GCM 
O'*' ground state might appear counter-intuitive, as, at 
small radii < 3 fm, the protons are obviously shifted to 
the inside. At larger radii, however, the tail of the den- 
sity of the GCM 0+ ground state becomes slightly larger 
than the density of the other states, which is almost un- 
detectable on the linear scale of Fig. [71 Because of the 
factor r"* in the mean-square radius integral in polar coor- 
dinates, this tail is much more important than the center 
of the nucleus. 

Figure. [7] puts into evidence that the reduction of the 
depletion factor at each stage of the calculation is partly 
due to the reduction of the maximum density pmax Obvi- 
ously, shell effects can reduce the density at some radii, 
but also enhance it at others. This indicates that the 
definition of the depletion factor ^ contains an ambi- 
guity concerning the reference density. An alternative 
definition of a depletion factor could be 

jp Psat,r Pcent^T /i i \ 

-Tsat.T = , iJ-lj 

Psat,T 

where Psa.t,T with t ^ p, n, t is now the saturation 
value of the proton, neutron, and total density. For '^^Si, 
we have p^at.p = (14/34) x 0.16 fm'^ = 0.066 fm-^, 
Psat.n = (20/34) X 0.16 fm-3 ^ 0.094 im-\ and p.^t^t = 
0.16 fm~^, respectively. Unlike Fmax, this alternative de- 
pletion factor i^sat can also be used to quantify central 
bumps in the density distribution. 

The results of Figs. [5] and [7] using these two depletion 
factors are summarized in Tables |T] and [Hi The value of 
Psat,p is systematically smaller than the values of 
and, as a consequence, the values of i^sat for protons are 
smaller than those of F^ax- For neutrons, the value of 
^max.ri IS always negative because of the central bump of 
the neutron density distribution. Again, the large central 
bump predicted by the HF calculation is reduced when 
correlations are added. Altogether, the central total den- 
sity is always larger than saturation density as evidenced 



TABLE II: Same as Table U but for the charge density, 
Eq. 



pccnt pmax -f^max -^'sat 

HF 0.056 0.071 0.21 0.15 

HFB+LN 0.060 0.071 0.16 0.09 

N&iZ 0.058 0.071 0.18 0.12 

N&iZ, J = 0.060 0.070 0.14 0.09 

GCM(g.s.) 0.064 0.070 0.09 0.04 



by Fsat,t. 

V. FURTHER DISCUSSION AND 
CONCLUSIONS 

There are two major differences between the structure 
of ■^■*Si as predicted by our calculation and the structure 
of other candidates for bubble structure that have been 
discussed in the 1970s-1990s. First, as discussed above, 
the depletion of the central density in '^'*Si appears for 
the proton density only, whereas the total density has 
an almost flat distribution throughout the bulk of the 
nucleus. Second, the level ordering of bubble nuclei is 
usually different from the one of regular nuclei, which is 
not the case for "^^Si. Take, for example, the hypothet- 
ical bubble-type configuration of the '^^Ar discussed in 
Ref. [1^ . There, the sequence of single-particle levels is 
altered such that the 2si/2 level is pushed above the 1^3/2 
level for both protons and neutrons. By contrast, for ^^Si 
only the relative distance of levels is changed such that 
the Z = 14 gap opens up, cf. Fig. [2] In the absence of 
a bubble structure of the total density, and of the rear- 
rangement of shells that is typical for bubble nuclei, the 
predicted anomaly of the proton density distribution of 
■^^Si appears to be an example of a central depression of 
the density, as observed also for many other nuclides [sj, 
rather than a nuclear bubble. 

Taking into account correlations reduces the central 
depletion of the proton density in "^^Si, as expected from 
earlier studies of other systems. Our main findings are: 

(i) A HFB+LN calculation overestimates proton pair- 
ing correlations in '^^Si. Particle-number projection 
of the spherical state constructed with HFB+LN 
reduces the pairing correlations, such that the den- 
sity profile almost goes back to the HF one. Clearly, 
a treatment of pairing correlations beyond the 
mean field with exact particle-number projection 
is needed in this case. 

(ii) Fluctuations in quadrupole degrees of freedom 
strongly even out the fluctuations in the density 
proflle; the central depression and the outer bump 
of the proton density are reduced, as is the central 
bump of the neutron density. 



8 



(iii) The central depletion of the density is less pro- 
nounced when looking at the experimentally ob- 
servable charge density instead of the point proton 
density. 

While all of the above findings can be expected to be 
generic on a qualitative level, the quantitative increase 
of the depletion factor when going from spherical HF 
to full projected GCM might depend on choices made 
for the effective interaction. In particular, in view of 
the somewhat too large mixing that we find between the 
two lowest 0+ GCM states, our calculation might slightly 
overestimate the role of shape fluctuations in the ground 
state. 

As discussed in the introduction, a central depletion 
of the proton density of ^''Si has been suggested as an 
explanation for the reduction by about 0.6 MeV of the 
spin-orbit splitting of the neutron 3/2~ and l/2~ lev- 
els inferred from transfer reactions [s^, ■ One has to 
be careful about such conclusion. First, the connection 
between the centroids of the spectral strength function 
of combined one-nucleon pick-up and removal reactions 
and the eigenvalues of the single-particle Hamiltonian is 
mo del- dependent [67| . and when looking at the domi- 
nating fragments only, the comparison is far from clear. 
Second, as the symmetry-restored GCM method corre- 
sponds to a superposition of many states obtained with 
different mean fields, there is no straightforward proce- 
dure that would allow for a statement about effective 
single-particle energies based on the density profile from 
our calculation. The only meaningful way to compare 
with the data would be to perform the same kind of cal- 
culation for '^^Si and '^^S, which at present, however, is 
out of our reach. 

When discussing the spin-orbit splitting of the neu- 
tron 2p levels in '^'^Si, there is an additional complication 



that goes even beyond these considerations. At spher- 
ical shape, the neutron 2^3/2 and 2pi/2 levels are far 
above the Fermi energy and outside of the energy interval 
shown in the Nilsson diagram of Fig. [5J In our spheri- 
cal HF calculation of '^"'Si, the former is weakly bound 
at —0.56 MeV, whereas the latter is even unbound at 
-1-0.83 MeV. In such a situation, the coupling to the con- 
tinuum has to be carefully taken into account, which is 
a task that goes beyond our study. By contrast, both 
single-particle levels are predicted to be (weakly) bound 
in a similar HF calculation for '^^S. 

In summary, we find that correlations from pairing and 
fiuctuations in quadrupole deformation substantially re- 
duce the central depletion of the proton density in '^'^Si. 
The extension of our method to the calculation of tran- 
sition densities in the laboratory frame as obscrvables in 
inelastic electron scattering is currently underway jsoj . 
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